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Abstract 



o. 

^SJ ' Following the notion that the Titius-Bode rule (TBR) may also be applicable to some extra-solar planetary 

_H ' systems, although this number could be relatively small, it is applied to 55 Cancri, which is a G-type main-sequence 

jrt . star currently known to host five planets. Following a concise computational process, we tentatively identify four 

l~^ ■ new hypothetical planetary positions given as 0.081, 0.41, 1.51 and 2.95 AU from the star The likelihood that 

\^ \ these positions are occupied by real existing planets is significantly enhanced for the positions of 1.51 and 2.95 AU 

in the view of previous simulations on planet formation and planetary orbital stability. For example, Raymond et 

al. (2008) argued that additional planets would be possible between 55 Cnc f and 55 Cnc d, which would include 

planets situated at 1.51 and 2.95 AU. If two additional planets are assumed to exist between 55 Cnc f and 55 Cnc d, 

the deduced domains of stability would be given as 1 .3-1 .6 and 2.2-3.3 AU. The possible planet near 1 .5 AU appears 

<~i , to be located at the outskirts of the stellar habitable zone, which is however notably affected by the stellar parameters 

Oh' as well as the adopted model of circumstellar habitability. We also compute the distance of the next possible outer 

Q \ planet in the 55 Cnc system, which if existing is predicted to be located between 10.9 and 12.2 AU, which is 

^ ■ consistent with orbital stability constraints. The inherent statistical significance of the TBR is evaluated following 

lyj \ the method by Lynch (2003). Yet it is up to future planetary search missions to verify or falsify the applicability of 

C^ i the TBR to the 55 Cnc system, and to attain information on additional planets, if existing. 

Key words: astrobiology — methods: statistical — stars: individual (55 Cnc) — stars: planetary systems 

QQ ' 1. Introduction tistical analyses. For example. Lynch (2003) argued that it 

f^ \ is not possible to conclude unequivocally that laws (or rules) 

f^ ■ One of the more controversial aspects in the study of the of Titius-Bode type are, or are not, significant, a conclusion 

l/~) ' Solar System as well as studies of selected extrasolar plane- that also appears to be consistent with the work by Hayes & 

r^^ . tary systems is the application of the Titius-Bode rule (TBR). Tremaine (1998). Moreover, the main conclusion of Lynch 

C^ ' Historically, the TBR was first derived and applied to the Solar (2003) is also shared by Neslusan (2004). Neslusan (2004) also 

^-H , System where it played a significant role in the search for new pointed out the curiosity that if the Earth's distance is excluded 

'""' ' planetary objects (e.g., Nieto 1972). The discovery of Uranus from the sequence, then the probability of the occurrence of 

^ \ by Herschel in 1781 and the largest object of the Mars-Jupiter that sequence by chance is radically reduced from 95 to 100% 

• '~j ■ asteroid belt, Ceres, by Piazzi in 1801 appeared to confirm the down to 3%. 

rN [ applicability and significance of the TBR. The TBR historically If no appropriate physical explanation for the TBR is identi- 
ty , given by the formula fied, it is typically argued that the TBR may be merely a rule of 
■■■' —O/i^nqo" — nioQ n\ chance ("numerology"), fuelled by the psychological tendency 

7*77 - — U.4I: ~\~ U.O X Z . Ill — - CX) ,U,-L,Z,0,... \^) •! • r 1 -/XT ^ i r\r\ a\ y~v 

to identity patterns where none exist (Newman et al. 1994). On 
(in AU) is, however, entirely inapplicable to Neptune and to the the other hand, there are previous detailed astrophysical and as- 
former Solar System planet Pluto', which is still often consid- trodynamical studies that potentially point to a physical back- 
ered to be an appropriate representative of the steadily increas- ground of the TBR, at least for the Solar System and (by impli- 
ing number of Kuiper Belt Objects (e.g., Jewitt & Luu 2001). cation) for a selected number of extrasolar planetary systems as 
A considerable amount of previous discussion revolves well. For example. White (1972) argued that jet streams may 
around the interpretation of the TBR for the Solar System, develop in a rotating gaseous disk at discrete orbital distances 
including its well-known insufficiency, see, e.g., Graner & given by a geometric progression. This result appears to be 
Dubrulle (1994), Hayes & Tremaine (1998), Lynch (2003), broadly consistent with the work by Schmitz (1984), who pur- 
and Neslusan (2004) for previous studies largely based on sta- sued an analytic analysis of the structure, stability, and form of 
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marginally stable axisymmetric perturbations of idealized ro- 

Pluto is mentioned merely for historical reasons; it is often considered to ... i j it- * j u a ^i. ^ ^u j- i » f tu 

,.,,., , ,. . . ■'_ , ^„„ i- , ^ , o .,, ■ ■ , tating gas clouds. His study showed that the radial parts of the 

highlight the limitations of the 1 BR tor the Solar System. Pluto is in the ° ° ■' ^ 

3:2 mean-motion resonance with Neptune and is thus not a dynamically expansion solutions obtained from the governing differential 

independent object. Moreover, long-term integradons have uncovered that yield simple oscillating functions with TBR-type features. 

its orbit is chaodc yet stable over biUion-year timescales (Malhotra 1993). A similar approach was undertaken by Graner & DubruUe 
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Table 1. Stellar Parameters 
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Quantity 



Value 



Reference 



Spectral Type 

Distance "^ 

Apparent Magnitude V 

RA Coordinate 

DEC Coordinate 

Effective Temperature 

Radius'^ 

Mass 

Age'' 

Metallicity [Fe/H] 



G8 V 

12.5±0.13pc 

5.96 

08'^ 52'" 35.8 IP 

+28° 19' 50.95" 

5280 K 

0.925 ± 0.023 i?0 

0.92 ± 0.05 A/q 

5 ± 3 Gyr 

+0.31 ±0.04 



Gonzalez (1998) 
ESA (1997) 
SIMBAD website^ 
SIMBAD website'' 
SIMBAD website*" 
see text 

Ribas et al. (2003) 
Valenti & Fischer (2005) 
Fischer et al. (2008) 
Valenti & Fischer (2005) 



"Data from the Hipparcos Catalogue. 

''See http : //simbad. u-strasbg. 



fr 



'^ Alternative value: 1.15 ± 0.035 Rq (Baines et al. 2008). 

''a more stringent determination, which is 4.5 ± 1 Gyr, has been given by Donahue (1998) and Henry et al. (2000). 

Table 2. Planetary System" 



Planet 


Distance (AU) 


Period (d) 


M sini (Mj) 


Discovery 


55 Cnc e 


0.01583 


± 0.00020 


0.7365400 ± 0.0000030 


0.027 


±0.0011 


2004/2011 


55 Cnc b 


0.115 


±0.0000011 


14.65162 ±0.0007 


0.824 


± 0.007 


1996 


55 Cnc c 


0.240 


± 0.000045 


44.3446 ± 0.007 


0.169 


± 0.008 


2002 


55 Cnc f 


0.781 


± 0.007 


260.00 ±1.1 


0.144 


±0.04 


2007 


55 Cnc d 


5.77 


±0.11 


5218 ±230 


3.835 


±0.08 


2002 



"Data given by Fischer et al. (2008) and references therein, except for 55 Cnc e, for which the data from Dawson & 
Fabrycky (2010) and Winn et al. (201 1) were used. 



(1994), who argued that a Titius-Bode type law emerges au- 
tomatically as a consequence of the scale invariance and rota- 
tional symmetry of the protoplanetary disc, and that such ge- 
ometrical relationships are a generic characteristic of a broad 
range of physical systems. Li et al. (1995) investigated the 
nonlinear development and evolution of density disturbances 
in nebular disks on the basis of one of their previous studies. 
Their analysis showed that the perturbed density is unstable 
with respect to self-modulation, leading to the formation of 
density field localization and collapse. In the case of a self- 
similar collapse, the perturbed density was found to increase 
with time and to form steady rings with very large but limited 
amplitudes. The distances of these rings were found to follow a 
geometric progression akin to the TBR. Also, in his summary 
about the final stages of planetesimal accumulation, Lissauer 
(1993) pointed out that "[t]he self-limiting nature of runaway 
growth strongly implies that protoplanets form at regular inter- 
vals in semimajor axis." This result appears to entail geometri- 
cal spacing of planets in planetary systems, potentially modi- 
fied by their long-term orbital dynamical evolution, in possible 
agreement with (a modified version of) the TBR. 

Another approach invoking statistical numerical experi- 
ments was adopted by Hayes & Tremaine (1998). They chose 
to fit randomly selected artificial planetary systems to Titius- 
Bode type laws considering a distance rule inspired by the Hill 
stability of adjacent planets. They did not identify a height- 



ened significance of the TBR, except that its meaning is that 
stable planetary systems tend to be spaced in a regular man- 
ner However, they indicate that their method could be used 
to identify potentially unstable planetary systems, especially if 
applied in conjunction with long-term orbit integrations. A fur- 
ther contribution to the interpretation of the TBR was made by 
Nottale (1996) and Nottale et al. (1997) based on an approach 
of scale relativity and quantization of the Solar System. It at- 
tempts to describe the Solar System in terms of fractal trajec- 
tories governed by a Schrodinger-like equation. The physical 
background and justification of this approach, however, remain 
highly uncertain. A controversial but testable aspect of this the- 
ory is that it proposes the existence of one or two small planets 
between the Sun and Mercury, which to date have no support 
through observations. On the other hand, these studies point to 
the general possibility of "empty" orbits, which may also be of 
interest to generalized applications of the TBR as well. 

In the following, we will apply a generalized version of the 
TBR, without subscribing to a specific physical interpretation, 
to the planetary system of 55 Cancri (= p^ Cnc = HD 75732 = 
HIP 43587 = HR 3522; G8 V) (see Table 1), which was his- 
torically the fourth star other than the Sun after 5 1 Peg, 70 Vir, 
and 47 UMa (not counting the pulsar PSR 1257+12 as well as 
other controversial cases) that was identified hosting a planet. 
A previous effort to apply the TBR to 55 Cnc has been pur- 
sued by Poveda & Lara (2008), which will also be discussed in 
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the following. Additionally, we will compare our results to a 
modified version of the TBR given by Lynch (2003) previously 
applied to the Solar System. 

55 Cnc is known to host five planets, named 55 Cnc b 
to 55 Cnc f, which were discovered between 1996 and 2007 
by Butler et al. (1997), Marcy et al. (2002), McArthur et al. 
(2004), and Fischer et al. (2008); note that all planets have been 
detected using the radial velocity method. The positions of the 
planets, according to their original determinations, are given as 
0.038, 0.1 15, 0.240, 0.781, and 5.77 AU, respectively, and their 
masses (A/pSini) range from 0.034 to 3.835 A/j (see Table 2; 
it also gives information on the respective uncertainties). 

However, a notable controversy emerged about the inner- 
most planet 55 Cnc e. Dawson & Fabrycky (2010) argued 
that the position of 55 Cnc e is given as 0.016 rather than 
0.038 AU based on an improved method that allowed to dis- 
tinguish an alias from the true planetary frequency in a more 
reliable manner This revised distance was subsequently con- 
firmed by Winn et al. (201 1) based on continuous photometric 
monitoring with the MOST (Microvariability & Oscillations of 
STars) space telescope. Moreover, there is an ongoing discus- 
sion about the principal possibility of additional planets in the 
55 Cnc system, particularly in the gap between 0.8 and 5 AU. 
This speculation is fuelled by detailed orbital stability studies, 
including studies involving putative Earth-mass planets (von 
Bloh et al. 2003; Rivera & Haghighipour2007; Raymond et al. 
2008; Smith & Lissauer 2009; Ji et al. 2009). If such a planet 
exists, it would also possibly imply the presence of a poten- 
tially habitable planet if appropriate conditions are met (e.g., 
Lammer et al. 2009). 

The star 55 Cnc is a middle-aged main-sequence star with a 
mass of approximately 0.95 Mq (Valenti & Fischer 2005), with 
a stellar effective temperature of about 5250 K (see Table 3). 
From the Hipparcos parallax of 79.8 ± 0.84 mas (ESA 1997), 
a luminosity of 0.61 ± 0.04 Lq (see Table 1) can be deduced, 
which is consistent with its spectral type G8 V (Gonzalez 
1998); see also discussion by Marcy et al. (2002). The age 
of the star can be obtained from the strength of the chromo- 
spheric Ca II Hh-K emission, indicating an age of 4.5 ± 1 Gyr 
(Donahue 1998; Henry et al. 2000); see also Baliunas et al. 
(1997) for previous work. Another element of our work is the 
assessment of 55 Cnc's circumstellar habitability, which will 
be pursued following the approach by Kasting et al. (1993) 
and Underwood et al. (2003); additionally, it will consider the 
possible extension of the outer boundary of habitability as dis- 
cussed by, e.g.. Forget & Pierrehumbert (1997), Mischna et al. 
(2000), and Halevy et al. (2009). 

Our paper is structured as follows. After the introduction of 
historical aspects of the TBR and 55 Cnc itself, as done, we de- 
scribe our methods and results. Emphasis will be placed on the 
calculational process concerning the TBR based on four differ- 
ent measures of difference; they are applied to assess the devia- 
tions between the positions of the five known planets of 55 Cnc 
(i.e., 55 Cnc b to 55 Cnc f) and the predictions of the TBR, thus 
allowing to constrain the unknown TBR parameters. In addi- 
tion, we comment on the possibility of a habitable planet in 
the 55 Cnc system as implied by the TBR. Thereafter, we con- 
sider other relevant studies, particularly studies regarding con- 
straints from orbital stabihty simulations and planet formation. 



(2) 



Table 3. Stellar Effective Temperature 

Tcff (K) Reference 

5150 ±75 Gonzalez (1998) 

5336 ± 90 Fuhrmann et al. (1998) 

5250 ±70 Gonzalez &Vanture( 1998) 

5243 ±93 deStrobeletal. (2001) 

5338 ± 53 Ribas et al. (2003) 

5279 ± 62 Santos et al. (2004) 

5234 ± 44 Fischer & Valenti (2005) 



We also comment on a previous appUcation of the TBR to the 
55 Cnc system given by Poveda & Lara (2008). Moreover, we 
provide an analysis of the statistical significance for the TBR 
forwarded by this work, as well as of the previous TBR ver- 
sion by Poveda & Lara (2008), following the method of Lynch 
(2003). Finally, we present our summary and conclusions. 

2. Methods and Results 

2.1. Computation of the TBR Parameters 

The centerpiece of our study is the application of the TBR to 
the planetary system of 55 Cnc. Akin to the application of the 
TBR to the Solar System, we originally selected as approach 

TB _ / j4 if i = 1 

"' ~ I A + B • Z'-2 if i > 1 

guided by the original planetary data of Fischer et al. (2008). 
Here A, B, and Z constitute free parameters, which need to be 
determined based on the known positions of the previously dis- 
covered five planets. Based on the five known planets 55 Cnc e, 
55 Cnc b, 55 Cnc c, 55 Cnc f, and 55 Cnc d, are identified as 
i = 1, 3, 4, 6, and 9, respectively. 

This version of the TBR is consistent with its original ver- 
sion customarily applied to the Solar System. Note that B in- 
dicates the compactness of the system (i.e., the higher B, the 
less dense the overall planetary distribution), whereas Z indi- 
cates the progression of the relative spacing of the planets. For 
the Solar System, the free parameters of the above given equa- 
tion take the following values: A = 0.4, B = 0.3, and Z = 2.0. 
However, in the meantime, a controversy emerged regarding 
the location of the innermost planet 55 Cnc e, which based on 
a re-analysis of existing data was identified to be positioned 
at 0.016 AU rather than 0.038 AU. We therefore revised the 
proposed version of the TBR such that the position of the in- 
nermost planet is changed to a™ = A/Z. Fortunately, test 
simulations have shown that the key predictions of this work, 
i.e., the location of the hypothetical planets at 0.081, 0.41, 1.51 
and 2.95 AU as well as the possible exterior planet, expected to 
be located between 10.9 and 12.2 AU (if existing), are entirely 
unchanged by this alteration. Moreover, the innermost planet 
is also not part of the polynomial fit determined by Z, just like 
in case of the Solar System. 

The pivotal aspect of our study is to determine and optimize 
the parameters A, B, and Z based on the positions of the five 
known planets of the 55 Cnc system, which are identified as 
i = 1, 3, 4, 6, and 9. The positions of these five planets are 
calculated using adequate trial values for A, B, and Z. These 
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Table 4. Titius-Bode Parameters 



MDIF 



A 



B 



Z 



MDIF 1 0.032 ± 0.003 

MDIF 2 0.031 ±0.001 

MDIF 3 0.031 ±0.002 

MDIF 4 0.031 ±0.001 



0.049 ± 0.001 
0.050 ± 0.001 
0.050 ± 0.001 
0.047 ± 0.001 



1.973 ±0.006 
1.971 ±0.006 
1.969 ±0.006 
1.991 ±0.004 



Table 5. Distances of Real and Hypothetical Planets" 



Planet 


Observed 




Titius-Bode Rule 








MDIFl 


MDIF 2 


MDIF 3 


MDIF 4 


55 Cnc e 


0.016 


0.016 


0.016 


0.016 


0.015 


... (HPL I) 




0.082 


0.081 


0.081 


0.078 


55 Cnc b 


0.115 


0.130 


0.129 


0.130 


0.125 


55 Cnc c 


0.240 


0.224 


0.224 


0.225 


0.218 


... (HPL 2) 




0.41 1 


0.411 


0.413 


0.403 


55 Cnc f 


0.781 


0.779 


0.781 


0.783 


0.771 


...(HPL 3) 




1.506 


1.509 


1.511 


1.504 


... (HPL 4) 




2.940 


2.945 


2.945 


2.965 


55 Cnc d 


5.77 


5.770 


5.774 


5.768 


5.872 



[Vol. , 



"Distances in units of AU. 




1 23456789 
Number of Planet 



Fig. 1. Distances of observed (full circles) and predicted (open circles) 
planets (or planetary positions) in the 55 Cnc system. The predicted 
planets have been computed following the TBR based on MDIF 1 (see 
Table 2). 



parameters are subsequently optimized through minimizing the 
numerical differences between the observed and "predicted" 
distances (i.e., semimajor axes) of the observed five planets. 
The later also requires a measure of difference. In the follow- 
ing four different measures of difference, henceforth referred 
to as MDIF, have been adopted, which are: absolute linear de- 
viation (MDIF 1), relative linear deviation (MDIF 2), absolute 
quadratic deviation (MDIF 3), and relative quadratic deviation 
(MDIF 4). The equations for the distance difference parameter 
5ri are thus given as 



„TB 



''p'.i '''pl,i 



rpi,i 



'pl,i 



5n 



Sn 



^pl.i 

( tb\ 

ypi.i — T^pi,ij 



^TB 



^pl,i 



(3) 
(4) 
(5) 
(6) 



where Tpj ^ denotes the observed distance (semi-major axis) of 
planet i and rZJ^ denotes the distance value obtained through 
the application of the TBR. 

Table 4 gives detailed information on the TBR parameters A, 
B, and Z, including the respective uncertainty bars, which are 
due to the distance uncertainties of 55 Cnc b to f (see Table 1) 
given by observations. It is found that A ~ 0.031, B ~ 0.049, 
and Z ~ 1.98. There is very little impact by the particular 
choice of measure of difference. For example, the spacing pa- 
rameter Z is found to vary between 1.969 ± 0.006 (MDIF 3) 
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12 3 4 5 6 7 8 9 10111213141516171819202122232425 

Serial Number of Simulation 

Fig. 2. Distance variations Sr for 55 Cnc f, 55 Cnc d, and the hypothet- 
ical planet 3 for a series of 25 statistical simulations. For 55 Cnc f and 
55 Cnc d, these variations are due to the observational uncertainties, 
which are 0.007 and 0.11 AU, respectively. For HPL 3, the distance 
variations were computed following the TBR based on MDIF 1 (see 
Table 2). 

and 1.991 ± 0.004 (MDIF 4). Furthermore, the distance uncer- 
tainties for 55 Cnc b to f, given by the observations, have also 
a relatively small effect on the values of A, B, and Z; note that 
the associated uncertainty is 0.3% or less in case of Z. 

2.2. Applications and Tests 

The derivation of the TBR parameters A, B, and Z can be 
used to identify possible additional planetary positions in the 
55 Cnc system (see Table 5), potentially occupied by plan- 
ets. Specifically, four planetary positions are found, which are: 
0.081, 0.41, 1.51 and 2.95 AU. As expected from our previ- 
ous discussion, the influence of the different choices of MDIF 
is only of very minor importance. Figure 1 offers a detailed 
depiction of the observed and predicted hypothetical planets. 
Note the very high quality of the fit for i > 1. 

Prior to discussing further results of our study, we want to 
present the findings from detailed tests about the convergence 
of the TBR parameters, notably Z. Again, the TBR parameters 
A, B, and Z are obtained by trying to fit the positions of the 
known planets 55 Cnc b to 55 Cnc f by using the TBR. This 
requires the choice of a measure of difference, MDIF; addi- 
tionally, it also requires the choice of a precision parameter for 
A, B, and Z that was chosen as lO^'^ with k — 1,2, 3, and 4. 
The control parameter of convergence was the summed resid- 
ual of the difference in the positions between the observed and 
predicted (fitted) distances of the five known 55 Cnc system 
planets. 

Our results are given in Table 6. As expected, it is found 
that the higher the selected precision, the more accurate are 
the obtained values for A, B (not shown), and Z deduced by 
fitting the positions of the five observed planets. However, the 
summed residuals are found to never approach exactly zero, 
regardless of MDIF and the precision parameter, although the 
residual will converge toward a very small limit (i.e., between 
10-2 and 5 x lO"'^), especially if MDIF 3 is selected. This 
result is expected. First, the higher the intended precision for 



A, B, and Z, the smaller the resulting residual. Second, the 
residual cannot approach zero, even if an increasingly higher 
precision for A, B, and Z is permitted, owing to the fact that 
the TBR (if applicable) should be viewed as a "rule" rather than 
a (physical) "law". There is no way that the TBR as proposed 
(or any modification thereof) will be able to precisely represent 
the positions of the currently known five planets in the 55 Cnc 
system as also encountered in applications of the TBR to the 
Solar System. 

Next we focus on the case of the hypothetical planet HPL 3 
(see Table 5) positioned near 1.51 AU; see Table 7 for de- 
tails. Concerning the four measures of difference, MDIF, the 
position of HPL 3 is found to range between 1.504 ± 0.025 
(MDIF 4) and 1.511 ± 0.028 (MDIF 3). Note again that the 
uncertainty in the position of HPL 3 due to the choice of MDIF 
is insignificant. However, this hypothetical planet is of poten- 
tially high relevance because it may be located in 55 Cnc's hab- 
itable zone (HZ). Therefore, it will be discussed in more detail 
below (see Sect. 2.3). 

Considering the potential importance of HPL 3, if existing, 
we also investigated the distance variation Sr due to the mea- 
surement uncertainties for 55 Cnc f and 55 Cnc d (see Table 1) 
as published; the uncertainty measurement is most substantial 
for 55 Cnc d (both in absolute and relative units). The results 
are given in Fig. 2; here the implemented version of the TBR 
was based on MDIF 1. Various positions for 55 Cnc f and 
55 Cnc d were assumed using a Monte-Carlo approach (Press 
et al. 1989) considering their observationally deduced distance 
measurements. We pursued a total of 25 simulations. In our 
series of simulations, the absolute variation 6r for HPL 3 was 
less than 0.03 AU. 

Another application of the TBR is the computation of the 
distance of the next possible outer planet of the 55 Cnc system, 
i.e., located beyond 55 Cnc d, which is at a distance of approx- 
imately 5.77 AU. This hypothetical planet is predicted to be lo- 
cated between 10.9 and 12.2 AU (see Table 8). More precisely, 
its position is predicted in the range between 11.33 ±0.43 
(MDIF 3) and 11.66 ± 0.39 (MDIF 4). 

2.3. The Habitable Zone of 55 Cnc 

In the following we discuss the possible location of a newly 
predicted planet in the HZ of 55 Cnc. The extent of 55 Cnc's 
HZ can be calculated following the formalism by Underwood 
et al. (2003) based on previous work by Kasting et al. (1993). 
Underwood et al. (2003) supplied a polynomial fit depending 
on the stellar luminosity and the stellar effective temperature 
that allows to calculate the extent of the conservative and the 
generalized HZ. Noting that 55 Cnc is less luminous than the 
Sun, it is expected that its HZ is less extended than the solar 
HZ, for which the limits of the generalized HZ have been given 
as 0.84 and 1.67 AU, respectively (Kasting et al. 1993). 

The luminosity of 55 Cnc has been deduced as 0.61 ± 0.04 
Lq based on the stellar distance and apparent magnitude; note 
that this value is fully consistent with the standard values for 
the stellar effective temperature and luminosity (see Table 1). 
However, if its stellar radius is 1.15 ± 0.035 Rq, obtained 
by Baines et al. (2008) based on new interferometric measure- 
ments, instead of R = 0.925 Rq (Ribas et al. 2003), the revised 
stellar luminosity is thus given as 0.90 Lq. This increased lu- 
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Table 7. Distance of the Hypothetical Planet 3 
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Fig. 3. Distances of 55 Cnc f (full circle) and the hypothetical planet 
3 (open circles) in relationship to the HZ of 55 Cnc, where dark gray, 
medium gray and hght gray denote the conservative, general, and ex- 
treme HZ, respectively (see Table 9). The position of HPL 3 was com- 
puted following the TBR based on MDIF 1, 2, 3, and 4 (from top to 
bottom). The depicted uncertainty bars ai'e 3cr, and are mostly due to 
the small observational uncertainties of 55 Cnc f and 55 Cnc d. The top 
figure shows the HZ of 55 Cnc based on R = 0.925 Rq (Ribas et al. 
2003), whereas the bottom figure assumes ij = 1.15 Rq (Baines et al. 
2008), implying a relatively high stellar luminosity of 0.90 Lq. 



Table 6. Convergence of Titius— Bode Parameters 



MDIF 


Precision 


Z" 


Residual 


MDIFl 


10-1 


2.1 


4.41 X 10-1 




10-2 


1.98 


3.90 X 10-2 




10-3 


1.973 


3.31 X 10-2 




10-" 


1.9704 


3.03 X 10-2 


MDIF 2 


10-1 


2.0 


3.93x10-1 




10-2 


1.98 


2.04 X 10-1 




10-3 


1.971 


1.91 X 10-1 




10-4 


1.9704 


1.90 X 10-1 


MDIF 3 


10-1 


2.1 


6.16x10-2 




10-2 


1.98 


5.94x10-4 




10-3 


1.969 


4.56 X 10-4 




10-4 


1.9691 


4.52 X 10-4 


MDIF 4 


10-1 


2.0 


3.80x10-2 




10-2 


1.99 


1.68x10-2 




10-3 


1.991 


1.68x10-2 




10-'' 


1.9908 


1.68x10-2 



°An equivalent increase in precision is also attained for 
the parameters A and B. 



MDIF Distance (AU) 

MDIFl 1.506 ±0.027 

MDIF 2 1.509 ±0.029 

MDIF 3 1.511 ±0.028 

MDIF 4 1.504 ±0.025 



minosity is, however, difficult, though not impossible, to rec- 
oncile with the other parameters of 55 Cnc, such as its apparent 
magnitude, spectral type and distance. 

Based on the standard values for the stellar parameters of 
55 Cnc, the limits of the conservative HZ are given as 0.75 
and 1.08 AU, whereas the limits of generalized HZ are given 
as 0.68 and 1.34 AU (see Fig. 3 and Table 9). Surely, larger 
limits are attained if the higher value for the stellar luminosity 
is adopted. The underlying definition of habitability is based 
on the assumption that liquid surface water is a prerequisite for 
life, a key concept that is also the basis of ongoing and future 
searches for extrasolar habitable planets (e.g., Catanzarite et al. 
2006; Cockell et al. 2009). The numerical evaluation of these 
limits is based on an Earth-type planet with a CO2/H2O/N2 
atmosphere. 

We point out that concerning the outer edge of habitability, 
even less conservative limits have been proposed in the mean- 
time (e.g.. Forget & Pierrehumbert 1997; Mischna et al. 2000). 
They are based on the assumption of relatively thick planetary 
CO2 atmospheres invoking strong backwarming that is further 
enhanced by the presence of CO2 crystals and clouds. These 
types of limits can be as large as 2.4 AU in case of the Sun; 
however, they depend on distinct properties of the planetary 
atmosphere, and are thus subject to ongoing studies and con- 
troversies (e.g., Halevy et al. 2009). Concerning 55 Cnc, the 
revised outer limit can be as large as 1.85 AU for the standard 
stellar parameters and up to 2.31 AU for the increased value 
of the stellar luminosity. These limits can be considered in 
the assessment of the hypothetical planet near 1.51 AU (see 
Table 7). Note that the position of this hypothetical planet is 
not significantly affected by the choice of MDIF. Based on the 
standard parameters for 55 Cnc, this planet would be located in 
the extremely extended HZ. If the increased value of the stellar 
luminosity of 0.90 Lq is assumed, it would be located in the 
generalized HZ. A depiction of the different scenarios is given 
in Fig. 3. 

Another important aspect regarding the circumstellar hab- 
itability of 55 Cnc concerns the planet 55 Cnc f. It is a hot 
Neptune with a mass of about 15% of Jupiter's mass; therefore, 
it is almost certainly unfit to host life. 55 Cnc f is located at or 
near the inner edge of the 55 Cnc's HZ, especially if the stan- 
dard value for the 55 Cnc's luminosity is adopted. However, 
55 Cnc f may be an appropriate object for hosting one or more 
habitable moons. Previous studies on habitable moons of extra- 
solar planets have been given by, e.g., Williams et al. (1997) 
and Donnison (2010). 
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Table 8. Distance of an Hypothetical Exterior Planet 



Table 9. Habitable Zone of 55 Cancri 



MDIF 


Distance (AU) 


MDIFl 


11.35 ±0.45 


MDIF 2 


11.35 ±0.44 


MDIF 3 


11.33 ±0.43 


MDIF 4 


11.66 ±0.39 



3. Consideration of Other Studies 

3.1. Constraints from Orbital Stability Simulations and 
Comments on Planetary Formation 

Any planet proposed to exist in a star-planet system by the 
means of the TBR or otherwise needs to pass the test of orbital 
stability. Previously, detailed simulations of orbital stability for 
the system of 55 Cnc have been given by, e.g., Raymond et al. 
(2008) taking into account the five known planets of the 55 Cnc 
system. Raymond et al. (2008) presented a highly detailed 
study with particular focus on the effects of mean motion reso- 
nances due to the planets 55 Cnc f and 55 Cnc d. They showed 
that additional planets could exist between these two known 
planets, including hypothetical planets located near 1.51 and 
2.95 AU. Specifically, if two additional planets are assumed to 
exist between 55 Cnc f and 55 Cnc d, the domains of stabil- 
ity are given as 1.3-1.6 and 2.2-3.3 AU. This is an important 
finding regarding the hypothetical, potentially habitable planet 
HPL 3, proposed at 1.51 AU (see Table 7). 

A subsequent study of orbital stability was given by Ji et 
al. (2009). They also considered the impact of mean motion 
resonances and identified various unstable locations. However, 
they identified a wide region of orbital stability between 1 .0 
and 2.3 AU, a potential homestead of habitable terrestrial plan- 
ets. These potential planets were also identified to have a rela- 
tively low orbital eccentricity. Furthermore, the orbital stability 
simulations by Raymond et al. (2008) and Ji et al. (2009) are 
also lending indirect support to the principal possibility of an 
exterior planet, which, if existing, should be located between 
10.9 and 12.2 AU from the star (see Table 8). 

Although orbital stability for the hypothetical planet near 
1.51 AU appears to be warranted, it is still important to gauge if 
such a planet could have formed in the first place noting that the 
existence of close-in giant planet(s) have been viewed to have 
an adverse effect (Armitage et al. 2003). Subsequent studies 
by Raymond et al. (2005, 2006) have pointed out that the for- 
mation and habitability of terrestrial planets in the presence of 
close-in giant planets is indeed possible. 

Raymond et al. (2006) gave detailed simulations also en- 
compassing the system of 55 Cnc. Assuming that the giant 
planets formed and migrated quickly, they found that terrestrial 
planets may be able to form from a second generation of plan- 
etesimals. For 55 Cnc, objects with masses up to 0.6 M^ and, 
in some cases, substantial water content are able to form; these 
objects are found in orbit in 55 Cnc's HZ. This type of result 
is also consistent with findings from the quantitative numerical 
program by Wetherill (1996), who demonstrated that habitable 
terrestrial planets can form for a wide range of main-sequence 
stars, noting that the number and distance distribution of those 



Description 



Distance (AU) 



Standard Alternative 



HZ-i (general) 


0.68 


0.84 


HZ-i (conservative) 


0.75 


0.93 


HZ-o (conservative) 


1.08 


1.35 


HZ-o (general) 


1.34 


1.66 


HZ-o (extreme) 


1.85 


2.31 



planets is relatively insensitive to stellar mass. This work has 
been expanded by Raymond et al. (2006) and applied to the 
special case of 55 Cnc. They found that the majority of the 
simulations resulted in the formation of terrestrial planets posi- 
tioned in 55 Cnc's HZ, including planets near 1.5 AU from the 
star. 

The existence of additional planets, including habitable 
planets, in the gap between 55 Cnc f and 55 Cnc d, i.e., between 
0.781 and 5.77 AU, also appears to be consistent with the exis- 
tence of a protoplanetary disk estimated to have a radial surface 
density profile such as (7{r) = r"^/"^ (Fischer et al. 2008). They 
also discussed the overall surface density of the disk, which is 
estimated to be considerably higher than in the solar nebula. 
The formation of terrestrial planets^ in the gap is furthermore 
indirectly implied by the relatively high content of heavy ele- 
ments as indicated by the enhanced metallicity of 55 Cnc itself 
given as [Fe/H] ^ +0.31 ± 0.04 (Valenti & Fischer 2005). 

3.2. Previous Application of the TBR to 55 Cnc 

A previous attempt to predict additional planets in the 
55 Cnc system based on the TBR has been given by Poveda & 
Lara (2008). At that time all currently known planets 55 Cnc b 
to 55 Cnc f were already discovered. Poveda & Lara (2008) 
assumed a "law" of the form of 



,TB 



a ■ e 



Xi 



(7) 



(in AU) with a and A as free parameters. Taking into account 
the planets 55 Cnc e, 55 Cnc b, 55 Cnc c, and 55 Cnc f, located 
between 0.038 (originally determined position for 55 Cnc e) 
and 0.781 AU, and represented by i = 1 to 4, a and A were de- 
duced as 0.0148 and 0.9781, respectively. This modified TBR 
is successful to also "predict" the outermost planet 55 Cnc d 
at 5.77 AU, counted as i = 6, which however was well-known 
when Poveda & Lara (2008) proposed the new law. 

Making 55 Cnc d also part of the fit results in a small ad- 
justment of the parameters a and A, which then acquire the 
values of 0.0142 and 0.9975, respectively. The main aspect of 
the work by Poveda & Lara (2008) is to predict a new planet 
for i = 5 and possibly also an exterior planet for i = 7 as well. 
The corresponding distances of these hypothetical planets are 
given as 2.08 and 15.3 AU; the corresponding orbital periods 
of the planets are approximately 3.1 and 62 years, respectively. 
None of these putative planets are expected to be located in the 



^ Following the study by Smith & Lissauer (2009), any additional planets in 
the 55 Cnc system are expected to have a relatively small mass to warrant 
the system's orbital stability. The low mass of any possible planet can also 
readily explain why so far they were able to escape detection. 
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55 Cnc's HZ. Moreover, the work by Poveda & Lara (2008) 
assumes the planet 55 Cnc e to be located at 0.038 AU rather 
than 0.01583 AU as identified by Winn et al. (2011); note that 
Poveda & Lara published their work before the updated result 
became available. Thus, it is unclear if or how the "law" by 
Poveda & Lara (2008) could be altered to accommodate this 
finding. A less-than-ideal approach could be to assign i = 0.11 
(or rounded to i = 0.1) to the innermost planet 55 Cnc e. 

4. Evaluation of Statistical Significance 

An important aspect of this study is to evaluate the statistical 
significance of the TBR as deduced. In this case, we largely 
follow the method previously applied by Lynch (2003) to both 
the Solar System and the Uranian satellite system. In the first 
step, the mean square deviation is calculated. For the TBR 
proposed in the present work, we find 

^it = ^ E [log(«^ - ^) - (l°g B + {2.-2)- log Z)]'(8) 

i 

where a, constitute the observed distances for the planets i = 3, 
4, 6, and 9. For A, B, and Z, we adopt the average values at- 
tained by the sets of runs pertaining to the four different mea- 
sures of difference, MDIF (see Table 4). Furthermore, we pur- 
sued a sufficently large number N of statistical model runs, 
where for each model j we obtain 



x] 



\Y,[\og{a,^A)-{\ogB + {i-2 + ky,)-\ogZ)]\9) 



with yi as a random number in the range [—1/2, +1/2] and k 
taken as fc = 2/3 (Lynch 2003). For the TBR as proposed, we 
find Xfit = 8.47 • 10^^; moreover, the average value of Xj is 
given as 1.56 • 10" ^ with N = 10^. 

Equivalent expressions hold if the method suggested by 
Lynch (2003) is applied to the version of TBR previously given 
by Poveda & Lara (2008). In this case, we find 

l.r-r. .. ..^^2 



Xfit 



4^ 



[log ai- (loga + iA)]' 



for i = 2, 3, 4, and 6, and 



x| 



-E 

4^ 

i 



log flj - (log a + (i + kyi)\)] ' 



(10) 



(11) 



For the TBR by Poveda & Lara (2008), we find Xfit = 9.67 • 
10~^ and the average value of Xj is given as 2.15- 10"^. 

Evidently, each individual value of Xj can either be Xj ^ Xfit 
or Xj > Xfit- The amount of Xj values with Xj < Xfit gives 
the probability that the sequence as calculated by the TBR and 
the observed sequence can occur by change. Based on this ap- 
proach, it is found that the newly proposed TBR is significant at 
a level of 9 1 . 8 %, whereas the TBR proposed by Poveda & Lara 
(2008) is significant at a level of 96.5%. Thus, the analysis of 
statistical significance of these two different versions of TBR 
for 55 Cnc shows that both of them are inherently highly signif- 
icant albeit the level of significance for the version previously 
given by Poveda & Lara (2008) is slightly higher Moreover, 
the TBR attained by Poveda & Lara (2008) is based on two free 
parameters instead of three free parameters for the description 



of the planetary sequence. However, the newly proposed TBR 
has the attractive feature that it is of virtually the same form 
as that of the Solar System; see Eq. (1). Here the parameter 
B describes the compactness of the system-', which is much 
higher than for the Solar System. Furthermore, the parameter 
Z describes the relative geometrical spacing of the system; it 
is virtually identical to that of the Solar System, i.e., Z = 2.0. 
Perhaps surprisingly, the TBR appears to be about as appli- 
cable to the 55 Cnc system as to the Solar System gauged by 
whether or not the degree of similarity between the predicted 
and observed distance sequences can occur by chance (see 
Introduction for data on the Solar System). Thus, the appli- 
cation of the TBR to 55 Cnc appears to be justified. However, 
considerable care is required to draw final conclusions because 
the fitting of five known planets assuming four gaps in the se- 
quence is much easier than fitting the continuous sequence of 
nine distances, i.e., with the mean asteroidal distance included 
and Pluto omitted; see Neslusan (2004). It will be the task of 
future planetary search missions to inquire if the planetary sys- 
tem of 55 Cnc follows (1) the TBR of this work, (2) the TBR 
by Poveda & Lara (2008), (3) another version of TBR-type se- 
quencing, or (4) no TBR-type sequencing after all, as found for 
the vast majority of the currently known planetary systems. 

5. Summary and Conclusions 

The goal of our study is to provide an application of the 
TBR to the 55 Cancri star-planet system. It is clearly un- 
derstood that the TBR does not have the same standing as the 
various well-established methods for the search and identifica- 
tion of extrasolar planets; see, e.g., Jones (2008) for a detailed 
overview. However, as pointed out in the Introduction, the TBR 
appears to be successful (within limits) for a large segment of 
the Solar System. Based on principal considerations, it would 
therefore be highly extremely unlikely if there was no other 
planetary system for which a TBR-type planet spacing would 
be realized, although the number of those systems may be rel- 
atively small. 

Thus, the predictions conveyed in this study follow the hy- 
pothesis that the TBR is applicable to the 55 Cnc planetary 
system. This approach entails the prediction of four planetary 
positions interior to the outermost planet 55 Cnc d, which are 
given as 0.081, 0.41, 1.51, and 2.95 AU. Whether these po- 
sitions are "empty", filled with asteroid / comet belts (like in 
one instant of the Solar System, where Ceres is customarily 
viewed as substitute for a planet) or occupied by a developed 
planet is beyond the scope of this work; in this case additional 
astrophysical requirements need to be met. 

Pertaining to this latter point, the strongest case for addi- 
tional planets around 55 Cnc is given for the possible planets 
near 1.51 and 2.95 AU, noting that they would be most com- 
patible with detailed studies of orbital stability (Raymond et al. 
2008; Ji et al. 2009) as well as previous studies about 55 Cnc's 
proposed protoplanetary disk (Fischer et al. 2008). For exam- 
ple, Raymond et al. (2008) argued that additional planets could 



The increased compactness of the planetary system of 55 Cnc compared 
to the Solar System is already implied by the known five planets, which 
inspired coining the term "packed planetary systems" hypothesis (e.g., 
Raymond et al. 2005). 
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exist between 55 Cnc f and 55 Cnc d, including planets located 
near 1.51 and 2.95 AU in the view of long-term orbital stability 
simulations. Specifically, if two additional planets are assumed 
to exist between 55 Cnc f and 55 Cnc d, the domains of stabil- 
ity were derived as 1.3-1.6 and 2.2-3.3 AU. These two ranges 
are fully compatible with the results of the present study. 

Moreover, the possible planet near 1 .5 1 AU is of particuar 
interest as it appears to be located at the outskirts of the stel- 
lar HZ. The existence of this planet is consistent with both or- 
bital stability simulations and studies of the formation of ter- 
restrial planets for 55 Cnc. For example, calculations of terres- 
trial planet formation by Raymond et al. (2006) showed that the 
majority of their simulations lead to the build-up of terrestrial 
planets in 55 Cnc's HZ, including planets near 1 .5 AU from the 
star Nonetheless, the final verdict of potential habitability of 
this possible planet (if existing) depends on a sizeable number 
of factors, including the definition of the zone of circumstel- 
lar habitability, considering that the outer boundary of the stel- 
lar HZ is notably impacted by a variety of astrobiological as- 
pects (e.g., Lammer et al. 2009) including the structure, density 
and compositions of the planetary atmosphere (e.g.. Forget & 
Pierrehumbert 1997; Mischnaet al. 2000; Halevy et al. 2009). 

The TBR as implemented follows closely the previous ver- 
sion as applied to the Solar System with A, B, and Z as free 
parameters, which were obtained through careful statistical fit- 
ting considering the five known planets of the 55 Cnc system. 
Here we obtain A ~ 0.031, B ~ 0.049, and Z ~ 1.98. The 
parameter B indicates the compactness of the planetary sys- 
tem, which is much higher than for the Solar System as also 
indicated by the conventional version of Solar System's TBR 
{B = 0.3). The parameter Z indicates the relative geometri- 
cal spacing of the planetary system; it is virtually identical to 
that of the Solar System (Z ~ 2.0). We also computed the or- 
bital distance of a the next possible outer planet of the 55 Cnc 
system, i.e., exterior to 55 Cnc d. Our study predicts that if 
existing it would be located between 10.9 and 12.2 AU from 
the star This distance is consistent with previous orbital stabil- 
ity analyses, which predict the position of a possible exterior 
planet as beyond 10 AU (Raymond et al. 2008). 

In summary, it is obviously up to future observational ef- 
forts to verify or falsify the existence of any possible additional 
planets. From a general perspective, these efforts will assist 
in gauging the applicability of the TBR outside of the Solar 
System, and will help to discriminate between different ver- 
sions of the TBR if found applicable to the 55 Cnc star-planet 
system. 

The author acknowledges valuable comments by an anony- 
mous referee, which also resulted in the addition of Section 4. 
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